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Abstract 

We consider the heat-kernel expansion of the massive Laplace operator on the 
three dimensional ball with Dirichlet boundary conditions. Using this example, we 
illustrate a very effective scheme for the calculation of an (in principle) arbitrary 
number of heat-kernel coefficients for the case where the basis functions are known. 
New results for the coefficients Bs , B§ are presented. 
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It is important to know the explicit form for the coefficients in the short-time expansion 
of the heat-kernel, K(t), for some Laplacian-like operator on a (^-dimensional manifold 
M.. In mathematics the interest stems for example from connections between the heat- 
equation and the Atiyah-Singer index theorem 0, whereas in physics the interest in this 
expansion lies for example in the domain of quantum field theory where it is commonly 
known as the (integrated) Schwinger-De Witt proper-time expansion @. 

If the manifold M. has a boundary dAi, the coefficients B n in the short time expansion 
have volume and boundary parts H. Thus 



oo 

d 



K(t)~(47rty-2 £ B k t k (1) 

fc=0,l/2,l,... 

with 

B k = [ dVb n + [ dS c n . (2) 

JM JdM 

For the volume part effective systematic schemes have been developed (see for example 
@). The calculation of c n is in general more difficult. Only relatively recently the coef- 
ficient C2 for Dirichlet and Neumann boundary conditions have been found || |6j. When 
using the general formalism of ref. f| for higher spin particles, Moss and Poletti |7j found 
a discrepancy with the direct calculations of D'Eath and Esposito || (see also 0). The 
latter results have been confirmed in ]nj where a new systematic scheme for the cal- 
culation of C2 has been developed in the context of the Hartle-Hawking wave-function 



of the universe for the case when the full set of basis functions is known [10]. Finally. 



very recently the discrepancy have been resolved [HJ and now the results found using the 



general algorithm |12| are in agreement with the direct calculations || |^, [T0|[ . 



We will use a variant of the approach of ref. JT0| in order to show that higher coefficients 



c n may be calculated very effectively for the case that the basis functions are known. To 
illustrate the method in this letter we will concentrate on the calculation of the heat-kernel 
coefficients of the elliptic operator 

(-A + m 2 )(f) ni i tm = (A^ /m + m 2 )(j) n j tm (3) 

in the three dimensional ball defined by Br = {x G 1R 3 , \x\ < R}. We will treat explicitly 
Dirichlet boundary conditions, <fi n ,i,m(\x\ — R) — 0. As will be clear afterwards, other 
boundary conditions and higher dimensional balls may be treated in exactly the same 
way. 

Starting point is the equation 

J l+ i(X n ^ m R) = (4) 

for the eigenvalues A„^ m , which are {21 + l)-times degenerated. We are especially inter- 
ested in the calculation of the heat-kernel coefficients defined in equation ([[]). Instead of 
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calculating the heat-kernel coefficients itself, we will concentrate on the zeta function of 
the operator, eq. and recover the heat-kernel coefficients using the relations [[0| 

Res C{s) = ( . B X~° : (5) 
(47r) 2 F(s) 



for s = f , z*=±, §; -^±1, for Z e JN , and 



C(-P) = (-WjfTEr (6) 
(47TJ 2 

for p G 1N . 

Using relation (f|) for the eigenvalues, one may write the zeta function in the form (for 
a similar procedure see Hi~4"| ) 



oo oo £ 

COO = EE E 

n=0 i=0 m=—l 

oo 



» 2 ' + 1 )/^(* 2 + "'V^ln7, +i (^), (7) 

where the contour 7 is counterclockwise enclosing all eigenvalues which are known to be 
situated on the positive real axis. As it stands, the representation (0) is valid for 3?s > |. 

Before considering in detail the /-summation, let us first construct an analytic contin- 
uation of the fc-integral in equation fl7|) alone and let us define for that reason 

dk , , o ov „ d 



with v — 1 + 1/2. Deforming the contour to the imaginary axis, the analytic continuation, 

valid in the strip ~ < 3?s < 1, may be found. A similar representation valid for m = 
has been given in [15]. 



In order to continue, the idea is to make use of the uniform expansion of the Bessel 
function I u (k) for v — > 00 as z = kjv fixed flEfl. Actually we use the expansion 



d 1 00 d (t) 

^-\nlJk)~ / (10) 
with £ = -^=s ■ Here the functions d n (t) fulfill the recurrence relation 



1 / d \ 1 n ~ 1 

d n (i) = - 1 2 ) - lj - - xjd fe (t)d n _ fe (t), (11) 
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starting with do(t) = 1. In order to calculate up to B§ one needs the first eleven coefficients 
d n , which can be easily calculated using the recurrence. 

Adding and subtracting the leading terms of the asymptotic expansion for v —>■ oo, 
eq. may be split into two pieces, 

N+l 

Us) = N u (s) + J2 KM, (12) 
i=i 

with 

^w-=p / * im 'l^w-i^) as) 




and 



raS/i/ \ v 7 / 

As it stands, the A z u , equation ([14]), are well defined (at least) in the strip 1/2 < 5Rs < 1. 
However, the analytic continuation in the parameter s to the whole complex plane in 
terms of known function may be provided. To explain afterwards some details of the 
calculation let us give explicitly only the A l v of the two leading terms in the asymptotics 

(0), 

Al(s) = ^^ Rm 

27T2 

x r(.-I)r(l-.),F I (-i,.-i : i;-( ; ^)'), (15) 
m~ 2s sin(7rs) 



47T 

v \ 2N 



Al(s) = 

x r(s)r(l - s) 2 Fi |l, s; 1; - (^^J J ■ (16) 

Similar expressions for higher AJ, can be calculated, e.g. using standard integration pack- 
ages. As mentioned, we did explicitly calculate the coefficients up to B 5 and thus needed 
A l u for 2 = 1, 11. Here 2-^1 (o> b, c; z) denotes the hypergeometric function 

The representation (O) has the following very important properties. First of all, by 
considering the asymptotics of the integrand in equation (|T3"D for k — > mRjv and k — ► 00, 
it may be seen that 

00 

iV( S )=£(2Z + l)iV; + l( S ) 

1=0 
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is analytic in the strip 1 — N/2 < 3?s < 1. For that reason it gives no contribution to the 
residue of £(s) in that strip. Furthermore, for s = —k, k = 0, 1, 2, 3, we have N(s) = and 
thus it does also not contribute to the values of the zeta function at those points. Together 
with eqs. (|5]), this yields, that the heat-kernel coefficients are only determined by the 
terms Ai(s) with 

oo 

A t (s) = Y,(2l + l)A] +l (s). (17) 

1=0 

However, the Ai(s) may be given in terms of Hurwitz zeta functions and an explicit 
representation of Ai(s) showing the meromorphic structure in the whole complex plane 
may be given. The sum in ( |T7D may be easily done by means of the Mellin-Barnes type 
integral representation of the hypergeometric function, 

r(c) 1 f, V(a + t)V(b + t)V(-t), , t 



2 F l (a,b;c;z) = v / — / dt -± L\ A ' " ' H0% (18) 

1 (ajl (b) 2m Jc 1 (c + £J 

where the contour C is such that the poles of T(a + 1) and F(b + t) lie to the left of it and 
the poles of r(— t) to the right ]T7|]. Defining 

h(a,b,c;n) ^f^(l + l) n 2 F 1 (a,b;c;- ) (19) 



1=0 



2/1 V A* 



and closing the contour C to the left, we arrive at 

oo (_-i)k 

h(a,b,c;n) = ^ i_L^* x (20) 

fc=0 fc! 

^ r ( ^. : ^ + , CH( _ n + 2a + 2 ^ 

which may be used for all summations we need for the calculation of the A^s. For example 
we obtain 

Ai = ™^M Bmr ( s -i)r(i- s ) ft (-i, s -i,i;i 

= R M r(1 _ f izl)! (mfl) ^£fci) c ( a + 2s _ 2; i 

2tt2 ^ Z! s + Z V 2 

* = - m ' 2 ' 9 sinW r W r(i- s )Mi, s ,i;i) 

27T 

= -^^^r(i - a) £ t^Wf r (z + 8 ) Ch (21 + 28-1-1 

2tt l=Q Z! V 2 
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and similar expressions for the other A^s, % = 1, 11. The sums appearing in the Ai(s) 
are convergent for \mR\ < 1/2. Using this representation, equations @, @, together 



with 

i\ i 



C//(l + e;-) = - + O(e ) 



r(e-n) = ^ + 

the heat-kernel coefficients may be easily determined. Summarizing, we find the following 
new results for the coefficients Bs , B 5 , 
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m 2 1 



4 D 2 



B 5 = 7T2 - m 4 i? 

2 V 6 120.R 2 

_ / 64 16m 2 4m 4 i? 2///7? :! 

3 71 V 9009i? 3 + 315# + 3~ 
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m 4 47 m 2 



5 5 n2 \ 12 20160i? 4 + 120i? 2 

/ 202816 64m 2 
B 4 = 7T -— — — — - + 




72747675# 5 9009.R 3 315i? 

4m 6 R m 8 R 3 \ 
+ 



18 J 

521 47m : 

Bl = 7T2 — + 




" I 36 443520i? 6 20160.R 4 

m 8 R 2 \ 



2A0R 2 12 J 

( 25426048 202816m 2 32m 4 

B 5 = 7T-— ————= + 



15058768725i? 7 72747675# 5 9009i? 3 

3m 6 m 8 i? m 10 i? 3 \ 
+ 



9A5R 9 90 J 

For m = the coefficient -B5/2 agrees with the one found by Kennedy [ |18|] . 

As mentioned, other boundary conditions, higher dimensional balls, and even higher 
coefficients may be found in exactly the same way without any additional complication. 
This and more details of the calculation will be given in a separate publication. 
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